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ABSTRACT 


Selected results of a three year project to study the climatological 
and micrometeorological conditions in a tropical evergreen forest, 
are summarized in this report. The forest area under study is in the 
interior of Thailand and is influenced by a cool, dry northeast monsoonal 
flow in the winter and a warm, moist southwest monsoonal flow in the 
summer. Information was collected during a three year field program 
sponsored by the U.S. Army Natick Laboratories and included ground 
stations temperature, relative humidity, precipitation, evaporation, 
hours of sunshine and tower measurements of temperature, dew point 
temperature, wind velocity, incoming and reflected solar radiation, 
incoming and outgoing infra-red radiation, spectral radiation and sub- 
surface temperature profiles in the forest and in a cleared area within 
the forest. 

Since the topics investigated during this three year period were 
of a diverse nature and too numerous to account for in detail in a single 
report, only a selected number of research topics will be presented here. 
The emphasis will be on the microclimate as evidenced by the temperature 
and velocity fields above and within the forest and clearing, their 
average seasonal characteristics, their variability, the stability 
conditions and the micrometeorological parameters of the site. 

To put the findings in a proper framework, the general climate at 
the site will be briefly outlined. The results are purposefully 
presented in an expanded format so as to enable the reader to easily 
extract useful information. 


A selected listing of a representative data base will appear in 


Appendix C. 


INTRODUCT' ON 


This report serves as a final report on Grant DRXRD-GS-13660 from 
the Army Research office, Durham, N.C. to the Department of Meteorology, 
University of Maryland. The objective of the research was to prepare 
an analysis of the climatological and micrometeorological data obtained 
during a scientific field program called TREND (Tropical Environmental 
Data). The TREND field experiment was under the sponsorship and manage~ 

ment of the Earth Sciences Laboratory, U.S. Army Natick Laboratories* 
and carried out by the Applied Scientific Research Corporation of 
Thailand in collaboration with several agencies of the Thai Government. 
The experiments were conducted during the years 1967-1970 as an inter- 
disciplinary study of a tropical dry evergreen forest environm it and 
embodied surveys of soil, vegetation and surveys of micrometeorological, 
hydrological, and biological conditions. 

The experimental design for the micrometeorological survey called 
for the collection and processing of temperature, dew point temperature, 
wind speed and direction at numerous levels along both towers, various 
radiation fluxes at the top of the towers and near the ground, and sub- 
surface temperature profiles near both towers. Six satellite ground 

stations recorded temperature, humidity, precipitation, evaporation, and 
sunshine duration. 

The results presented in this report emanate from efforts at the 
University of Maryland to extract useful climatological and micrometeoro- 


logical information from the raw data available to us. 


*(Dr. Paul Dalrymple: Project Supervisor). 


l. Site and instrumentation 


The study site for the TREND experiment is the Sakaerat Forest, 


an area cf 


about 80 square kilometers in Thailand situated on the Korat 
Plateau at approximately 14931'N, 101955'E which is about 60 km south- 
southwest of Khorat (Nakhon Ratchasima) and about 190 km northeast of 
Bangkok (Figure 1). Meteorological data were taken along two towers 
approximately 50 m high (Fig. 2) and 450 m apart, one tower in the 
forest, the other in a small clearing and both based at about 535 m 
above sea level on a tilted slightly dissected sandstone plateau which 
slopes from southwest to northeast. Figure 3 shows schematically the 
location of the various measuring instruments along each tower. The 
forest includes good stands of dry evergreen forest, dry dipterocarp 
forest and clearings in various stages of regeneration. The immediate 
vicinity of the towers is comprised of a two-storied forest. The top 
story consists almost exclusively of Hopea ferrea, with a scattering of 
Shorea sericeiflora, with dense and continuous crown canopy with tree 
tops ranging from 20-35 m high. The second story consists of the species 
Hydrocarpus ilicifolius, Walsura trichostemon, Aglaia pirifera, 
Lagerstroemia, and Memecylon ovatum with tops ranging from 5-17 m high. 
The forest floor was well covered with thick undergrowth. The adjacent 
experimental clearing site was approximately 500 m in diameter, covered 
with low vegetation (mainly elephant grass), at different stages of 


regeneration. 


A general account of the instrumentation and monitoring procedures 
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is given in ASRCT (1969) and Dalrymple (1975). Here, only a brief outline 
of the instruments and data acquisition system will be presented. 
Temperature measurements at both towers were made with Hewlett - 
Packard quartz thermometers. The sensing element is a quartz crystal 
whose temperature dependent probes are placed in shielded, motored 
aspirators supplied by Climet (Fig. 4)- Two types of temperature probes 
were used, one for air and one for soil. 
Dewpoint temperatures were measured by a dewpoint hygrometer composed 
of two units: the mirrored transducer and the controller. A light 
beam was sensed by two separate photodiodes; one sensed the direct 
beam from the lamp, and the other sensed the reflected beam off the 


mirror. The two energy levels were continuously compared and when 


dew formed, the transducer output was applied to the controller circuitry 
where it was compared to a fixed voltage. Windspeed sensors were 

Climet three-up anemometers (Fig. 4) having a threshold velocity of 

0.3 m/sec. Wind direction was detected by Climet wind direction trans- 
mitters and vanes. 

Rainfall was measured electronically by a standard weighting-type 
gauge from Science Associates. The gauges were located at the top and 
near the base of each tower. 

Atmospheric pressure was measured with a precision pressure gauge 
made by Texas Instruments. 

There was in operation an extensive solar radiation program which 
included arrays of instrumentation mounted on the top of the forest 4 


tower and on the forest floor. There were also instruments placed 
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near the surface close to the tower in the clearing. Global radiation 
was measured by Eppley precision spectral pyranometers. Additional 
pyranometers measured the standard WMO broad band spectral regions with 
Schott colored glass filters: 
OG1 - 280 - 2800 nm 
RG2 - 630 - 2800 nm 
RG8 - 700 - 2800 nm 
GG14 - 475 - 2800 nm 
WG7 - 280 - 2800 nm 
Narrow band spectral regions were measured by interference type 
filters, designed to isolate relatively small wavelength intervals, 100 nm 
in the visible and several hundred nm in the near IR. The limits of 


the narrow band filters for the two TREND sets were (in nm): 


Above Forest Canopy Below Forest Canopy 
290-425 290-440 
335-505 380-540 
505-695 495-695 
650-1200 700-1050 
800-2350 800-2100 

1800-3000 1800-3000 


Complete sets of spectral pyranometers were mounted on the forest 
tower platform as well as below on the forest floor. Three normal incidence 
pyranometers measured the total direct short wave radiation as well as 


the spectral regions determined by the Schott OG] and RG2 broad band filters. 


Daylight illumination and UV radiation were measured with Eppley. 


photometers. Both utilized Weston selenium barrier-layer photocells. 
The wavelength response for the UV meter was restricted to 290-380 nm. 
Incoming and outgoing long-wave radiation above and below the canopy were 


measured by a Funk polyethylene radiometer. 
2. The data and their availability , 


The preparation of the raw data for computer processing was a 

formidable task because of the condition of the tapes. The major problem 
in reading the tapes was related to the incremental recorder's tendency 
to stretch the tape (from 2400 ft to 2700 ft ), resulting in poor spacing 
between characters. This spacing is interpreted as an inter-record gap, 
resulting in loss of data or loss of position on the tape drive. 
An older model tape drive (IBM 1401) was found to be able to read the 
tape better than an advanced tape drive. Subsequently, the data were 
repacked on a permanent storage tape in a compressed format compatible 
with Univac 1108 and CDC 6600, 7600 hardware. 

The technical details about the repacking and use of the repacked 
"D" and "A" tapes will be presented in a separate ‘User's Manual" which 
will outline the extensive development of data reduction schemes which 


were necessary to convert the raw data into a data resource useful for 


analysis. 


3. The climatic regime 


3.1. The regional climate 
Ramage (1971) defines a weusobdal region to be one in which: 

1. The prevailing re direction shifts by at least 120° 
between January and July. 

2. The average frequency of prevailing wind directions in 
January and July exceeds 40%. 

3. The mean resultant winds in at least one of the months 
exceed 3 m/sec. 

4. Fewer than one cyclone-anticyclone alternation occurs 
every two years in either month in a 5° latitude-longitude 
rectangle. 

Consequently, the monsoonal region would be bounded by the latitude- 
longitude rectangle 35°N - 25°s, 30°W - 170°E, and thus contain Thailand. 

The beginning of the monsoon season in inland Thailand is not 
clearly defined. In the peninsula of Thailand, Crutcher et al. (1969) 
define the onset and retreat of the monsoon as a period of 3 - 5 con- 
secutive days of rainfall (about May 5 and October 19). Usually, the 
Southwest monsoon is established in May and ends in September. The 
Northeast monsoon starts in November, and ends in February, with two 
transition periods. This pattern was also verified at the experimental 
site (Fig. 5), based on three years (1967-1970) of surface observations 
of climatological conditions. Accordingly, the selected periods for 


study span different stages of the monsson cycle and are of approximately 


two week time intervals. 


3.2 Climatological conditions at the experimental site 

To place the micrometeorological studies of the forest into a proper 
climatological context, data from the principal *surface station were 
summarized over the nearly three years of record to construct a climatology 
of the area. Computer printouts of hourly values of temperature, relative 
humidity, precipitation, pan evaporation and hours of sunshine were 
available from this field experiment. Also, derived climatological 
parameters were computed. 

Using the average air temperature and relative humidity values for 
each month, values of “effective temperature" were obtained in 
conventional fashion (Landsberg, 1970). Basically, the effective temperature 
takes into account the combined effect of temperature and humidity 
on a person's comfort. Effective temperature is defined as a sensation 
in which, for a given air temperature and humidity, the state of comfort 
is equal to that experienced for an environment at a lower temperature 
with saturated conditions. 

Precipitable water, U, in the atmospheric column over Sakaerat 
Experiment Station was computed using the method of Smith (1966) with 
empirical (A) values taken from Lettau (1970). These were computed 
using monthly average values of temperature and relative humidity. 

Values of the storage of water vapor in the atmospheric column were 
computed by taking the differences of precipitable water from month to 


month, AU. The advection, A, of water vapor into the atmospheric 


*There were 8 climatological surface stations in operation during the 


experiment (Fig. 6). 


column was computed as a residual in the balance equation, 


Storage = Evaporation - Precipitation + Advection 


AU = E - P +A 


Table 1 gives the three year average of climatological elements 
and the derived climatological elements for each aontR of record at the 
forest site. Figure 5 shows a representative annual variation of each 
element. The number of hours of daylight were taken from the 
Smithsonian Meteorological Tables on Duration of Daylight for the 15th 
of each month. 

During the cool, dry northeast monsoon, associated with high pressure 
over the Asian continent, the temperatures and relative humidities 
are comparatively low. The representative annual cycles of temperature 
and effective temperature show maximum values (27.5°C, 24.9°C) during 
the transitional month of March and minimum values (21.5°C, 20.0°C) 
in December. The annual average temperature at Sakaerat was 25.0°C 
and annual average effective temperature was 23.5°C, with only a 5-6C° 
variation in either during the year. The relative humidity generally 
falls during the winter monsoon period. The relative humidity is 
generally largest in October, with an annual average value of 85%, 
and smallest in February with annual average value of about 64%. The 
largest difference between air temperature and effective temperature 
occurred during February, 1969 when the relative humidity reached its 


period low value. The average air temperature during that month was 


27.2°C while the effective temperature in the relatively cool dry 
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monsoonal air was only 23.4°C, almost 4C° cooler. The northeast dry 
Monsoon occurs during a period of relative minimum daylength yet the 
sunshine recorders measured relative maximum amounts of sunshine during 
this period which indicates smaller amounts of cloud cover. The 
precipitation traces are also small averaging generally less than 1 mm/day 
during the northeast monsoon. 

The warm and moist southwest monsoon generally occurs during the 
period April to September corresponding to the months of maximum daylight. 
The temperature, precipitation and relative humidity are generally high 
during this period. The total rainfall during these six months was 
980 mm for 1969, 939 mm for 1970 and a relatively smaller amount of 
669 mm during 1968. The double maxima in precipitation during the 
southwest monsoon, shown in Fig. 5, is also evident in a summary 
of Bangkok climatological data for the period 1943-1952 from the Royal 
Thai Navy which was prepared by Lettau (1970). 

The relative humidities are high during these monsoonal periods yet 
the effective temperatures are lower than the air temperatures by a 
greater amount than during the dry northeast monsoon where the air temp- 
eratures are much lower. The humidity generally increases during the 
southeast monsoon period achieving its maximum toward the end of the 


period. 


The higher degree of cloudiness associated with the southwest 
monsoon is reflected in the sunshine recordings where relatively smaller 
amounts of sunshine occur each day than during the northeast monsoonal 


period. In fact, the monsoon cycle is well illustrated by the sunshine 


recording which is roughly 180 out of phase with the cycle of daylight 
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duration. 


The precipitable water record tends to parallel the relative humidity 
record over the whole period. The evaporation tends to respond to the 
solar radiation available at the surface and is also controlled by the 
ambient relative humidity and air temperature. The evaporation is 
generally largest when air temperatures are high and ambient humidity 
is low, and generally smallest during periods of higher ambient humidity, 
falling temperatures and fewer hours of sunshine. 

The water vapor storage term is generally small throughout the 
period with values within + 10 mm/month. The advection term shows the 
large monsoonal influx of moist air feeding the rainfall systems of 
May and September with outfluxes of moist air from the region during 
periods of evaporation excess over precipitation during the northeast 
Monsoon. Average rainfall during May and September is 172 mm/month 
and 252 mm/month respectively separated by three months of lesser rainfall 
with local minimum during July of 78 mm/month. The precipitation drops 
to near zero during the northeast monsoon season of December, January 
and February. For the rainfall peak of May, about 90% (146 mm) is 
contributed by evaporation, 17% (29 mm) is contributed by advection 
with 3.3 mm added to storage. During September, 62% (156 mm) is contri~ 
buted by advection, 37% (93 mm) is contributed by evaporation and 1.2% 
(3.0 mm) is contributed by storage. During the relative minimum of 
July, there is a fairly vigorous advection of about 43% (57 mm) of the 


evaporated water out of the region leaving only about 78 mm to rain out. 
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During the dry months of December, January and February nearly all the 
evaporated water is advected out of the column with little precipitation 
and only minor changes in storage. 

A more detailed description of the site, instrumentation and 


climate elements can be found in Thompson and Landsberg (1975). 


4. The microclimate of the experimental site 


The inventory of well established forest microclimates is 
meager. Relatively, little is known about the microclimates of tropical 
forests*. Perhaps their remoteness from research centers, and their 
adverse effect on the long term maintenance of the instrumentation, 
have restricted experimentation in such environments. The microclimates 
of clearings are also of interest. As already noted by Lee (1978), 
forest openings play important role in stand management and regeneration. 
The seedlings in a forest clearing are subjected to both higher 
radiation energy fluxes and surface temperature during the day and 
greater energy losses and lower temperatures at night. 

In this section, the microclimatological characterization of this 
dry (1500 mm annual zone) tropical forest and a nearby clearing is 
presented. This type of characterization supplements existing studies 
by virtue of its geographical location vegetation type, rainfall pattern 


and illustration of seasonal variability. 
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* Updated information and extensive bibliography on coniferous forests 
can be found in Jarvis and Landsberg (1976); on deciduous forests in 
Rauner (1976); Leigh (1975) addresses himself to the ecological aspects 


of tropical rain forests all over the world. 


4.1 Temperature time-height cross-sections in the forest 


The predominantly clear winters (8 hrs of sunshine) and cloudy 
summers (5 hrs of sunshine) typical of a monsoonal climate regulate 
to a large extent the temperature field at this Thailand forest. 

The two week averaged temperature time-height cross sections 
along the forest tower (Figs. 7, 8, 9, 10) suggest a subdivision of 


the canopy into three responding sublayers: 


(a) surface layer - 0 - 5m 
(b) sub~canopy layer - 5 -25m 
(c) canopy layer - 25 - 35m 


In (a) the influence of the ground is dominant, and as such, it 
is a region of frequent and strong stratified conditions. Ground 
moisture and cloud cover tend to subdue the extremities of these conditions, 
the length of their duration, and the depth of the ground affected 
layer. It is most pronounced during the dry and clear winter season. 

During July and September an isothermal structure prevails in 
sublayer (b). In January and April, isothermal structure forms primarily 
during the noon hours (10:00 a.m. - 4:00 p.m.). (Throughout this paper 
L.S.T. is used). During the remaining part of the day in January and 
April the (b) layer is characterized by inversions of 1-2°C/20 m. 

When clear conditions prevail, the canopy layer (c) could be viewed 
as the response layer to the solar heating. The intensity of this 
response differs from month to month, depending on additional ambient 


variables. For January, the temperature distortion in the top canopy 
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layer is strongest. This is a dry period with low winds (2.4 m/sec). 
The moisture stress and poor mixing might be responsible for this 
Pronounced distortion. In April, which is still relatively clear, 

the ambient winds increase (4.1 m/sec) and a such, the intensity of the 
upper disturbed layer decreases. During September (the September data 
were averaged over seven days only, resulting in less smooth cross 
sections as compared to others), the “above” canopy region (from 40 m 
up) is generally isothermal. The average temperature difference between 
the ground and canopy top around noon time is approximately 2°c in 

July and September, reaching 4°C in January and April. 

The previously discussed cross-sections attest to the dense nature 
of the canopy. Direct sunlight does not generally reach the ground 
(this was also verified from radiation measurements taken at the forest 
floor (Pinker et al. (1980))). The clear conditions of the dry 
winter monsoon are responsible for the much higher amplitude of the 
diurnal temperature wave at forest air interface. During the wet summer, 
the heating intensity is reduced at interface, and more homogeneous 


conditions prevail throughout the canopy. 


4.2 Temperature time-height cross sections in the clearing 


The structure of the temperature time-height cross sections in 
the clearing also calls for a three layer response subdivision: 

(a) surface layer - 0 -10m 

(b) “transition” layer - 10 - 25m 

(c) canopy affected layer - 25 - 35m 


The degree of distinction or even existence of these layers is 
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again a function of solar radiation intensity and mixing effectiveness. 

In January, the maximal temperature difference between ground and 
above canopy levels amounts to 2°c. From early evening hours till 
sunrise, inversions in the surface layer of 1°C/15 m build up and 
spread to the second layer (b). They are destroyed at sunrise and 
a lapse condition starts to build up, reaching 2°c/30 m at noon time. 
The existence of grass in the clearing up to a height of 1 m caused 
the displacement of the warmest layer above the ground (Fig. 7b, 8b) 
(it was subsequently cleared). During January, the edge effect of the 
forest canopy on the clearing is most pronounced, affecting a 10 m 
layer centered around the canopy top (at 30 m). 

In April, the difference between the maximum ground surface and 
above canopy temperatures is only 1°C due to the more turbulent structure 
of the wind field. Starting at the 2 m level, late evening inversions 
become established and expand in depth with time. After 8:00 a.m. 
lapse conditions prevail and increase in depth till noon time. The 
edge effect of the forest canopy diminished as compared to January. 
The "above canopy" region (from 40 m up) is generally isothermal. 

In July, the influence of the surface is confined to a shallow 
layer. The ground inversions are almost non-existent, while lapse 
conditions do exist in the first 4 m from 8:00 a.m. till 5:00 p.m. 
During the rest of the time, isothermal conditions dominate. This is 
a result of the more moderate mode of heating during the prevailing 
cloudy conditions and smaller radiative cooling during the night. 


In September, the lapse conditions close to ground are even milder 


= 
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than in July. The increased moisture content of the ground regulates 


the heating at the ground. 


4.3 Wind speed time~height cross sections in the forest and clearing 


A peculiarity in the surface winds over the Khorat site was noted 
(based on data for seven different months). During February, July 
(Figs. 11, 12) April, June, and September (not illustrated here) the 
maximum in wind speed occurs before sunrise and mid-night. The maximum 
in the surface boundary layer is generally observed at noon-time, as was 
the case in January and December (Figs. 13,14) at the Thailand site. 

A satisfactory explanation of this pattern has not yet been found. 

The wind speed cross sections along the forest and clearing towers 
will be illustrated for January, April and July (Figs. 15, 16, 17); 
for the other months the pattern repeats. 

A comparison of the forest and clearing average wind flows in- 
dicates that: 

a) the canopy strongly reduces the wind speed (for quantitative 

estimates see Table 2) 

b) the amplitude of the daily cycle is diminished over the forest 

as compared to the clearing 

c) during January, the amplitude of the daily cycle is enhanced 

close to the surface in the clearing (up to 16 m), while in 
April and July it is drastically suppressed below the 16 m 
level. (The diurnal heating cycle enhances the turbulent mixing 
in the clearing by day and results in higher wind speeds. In 


January, it is in phase with the above canopy flow and so amplifies 
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the effect. During periods with an above-canopy nocturnal 
maximum the heating effect is out of phase and dampens the 
cycle). 

d) in the forest, the amplitude of the daily cycle becomes sup- 

pressed below the 32 m level; the January case (Fig. 13) reveals 
a small increase in the wind speed around noon-time at the 
25 m level. 

According the Sadeh (1971), the one story forest canopy flow can 
be divided into three distinct zones: 

(a) flow within the trunk 

(b) flow within the tree crown - canopy flow 

(c) flow above the canopy. 

In this study, due to the two story nature of the forest, the 
distinction between (a) and (b) is not obvious. Seemingly a narrow 
"interaction" layer exists close to the canopy top and a "canopy flow" 
throughout the rest of the forested area. 

Figures 15, 16, 17 imply that the layer below the 30 m level is 
almost impermeable to the wind flow. Yet, there is a noticeable dif- 
ference between the three cases. The depth to which the | m/sec isotach 
penetrates does depend on the magnitude of the ambient wind speed. 

In a previous study, Thompson and Pinker (1975) computed the 
roughness parameter (2,) displacement height (d) and friction ve- 
locity (u,) for this site. They used a set of half hourly averaged 
wind speed profiles in the 30 - 46 m layer above the forest canopy 


which were classified as neutral. Table 3 presents a summary of the 
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roughness parameters, datum devel displacements and friction velocities 
as previously obtained. It seems that the values of 29 m for d in 
January and 27 m in July are consistent with the different depth of 
penetration of wind during those two periods. These rather small 
differences seem to have a significant effect on the diffusion processes 
in the canopy. The perturbed temperature structure at the canopy/air 
interface during January, is partially due to the low winds and poor 


mixing ( but not entirely, as the April case indicates). 
4.4 Subsurface temperatures in the forest and clearing 


Exchange of radiation at the ground causes periodic variations 
of the surface temperature. These changes affect the temperature of 
the soil. Flux of heat into and out of soil is a process of conduction 
and is proportional to the rate of change of temperature with depth. 
The constant of proportionality is the thermal conductivity and de- 
pends on such soil properties like porosity, moisture content and 
organic matter content of soil. The ground temperature will depend 
not only on the amount of heat transported to it but also on its 
ability to absorb that heat (thermal diffusivity). 

A detailed soil survey of thirty square kilometers of the forest 
experimental site was made by the Department of Land pavelopuent; Min- 
istry of National Development, Thailand. Soils were studied in borings 


and profile pits. Majority of the soils belong to the Red-Yellow 
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Podzolic great soil group; the parent material is sandstone (sand, silt, 
clay). Soil conditions in the core area are estimated to be fairly 
representative for general soil condition in the hilly parts of South- 
east Asia and those of Thailand in particular (for details see Bos (1968)). 
In all four cases studied, the amplitude of soil temperature 
decreases with depth in the bare and vegetation covered soil and the 
diumal wave penetrates deeper in the clearing than in the forest. 
At 40 cm the wave is almost damped out everywhere. In September due 
to the large water content of the ground the insulating effect is 
reduced, the heat conductivity is increased and the diurnal wave is 
found to penetrate deeper in the clearing. In the forest on the other 
hand, the effect of the increased heat capacity is dominant resulting 
in an almost isothermal structure of the soil. 
Generally, in dense stands, the diurnal and seasonal temperature 
waves are greatly damped in comparison to clearing and the heat flux 


to the ground is reduced. 


5. Microclimatic variability 


Information about the fluctuating nature of the temperature and 
velocity fields in the surface layer, over a rough terrain, is of 
significance to problems of local atmospheric diffusion, and for the 
assessment of the contributing factors to the generation of turbulence. 


Therefore, the available data, collected at the TREND site were analyzed 


Te te ee ene ee 
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to provide a better understanding of: 
a) the wind velocity and temperature spectra in the frequency 
range of 1072 ~ 102 cy/hr 
b) the dependence of turbulence spectra on stability in the 
frequency range 10 - 2 x 102 cy/hr 
c) the time-height dependence of the temperature and wind 
velocity variance 
In this report, only the time-height dependence of the temperature 
and velocity variances will be discussed, while a) and b) will be presented 
independently. 

Figures 19-22 represent the standard deviations (0) of temperature 
for January, February, April and July. They were obtained by computing 
half hourly o's from three hour averages over a period of two weeks. 
The o's for wind speeds were obtained in a similar manner and are 
presented in Figs. 23-26. 

During all the periods that were investigated, the largest o's 
of temperature were obtained before noon, during the steep change in 
temperature; rather than at the early afternoon when the temperature 
maximum is obtained. 

Inspection of the temperature cross sections in Figs. 7-10 attests 


to the following differences in the average maximum temperatures: 


Jan. Feb. Apr. July 
Forest (°C) 25 30 28 27 
Clearing (°C) 26 32 29 28 


Sects SSE ran ae ne ee ape nn SN aes Sie ea sat ere 
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The higher values of temperature and heating rate (clear skies) 
during February and April are probably the cause for the higher o's 
for this period. Figures 11-14 provide the basis for the differences 
between these two months. In February, from midnight till early morning, 
the winds are higher than in April. The mixing is more effective and 


the T's are lower. During the daytime, the February winds speeds 


drop below the April values, resulting in less mixing and higher o's. 

The July maximum temperature is by 2°C higher than the corresponding 
January value. Yet, the January Temperatures are higher, which might 

be -due to the lower wind speeds and the prevailing clear sky conditions, 
which create larger temperature gradients. 

The largest o's for the wind velocity occur during April and July 
when the winds are strongest. The maximum April value well fits the 
general pattern of the wind structure for this transition period in 
this monsoonal region. The lower July values agree well with the more 
persistent structure of the summer monsoon flow. 

In February, the largest o's occur closer to noontime. Fig. 11 
shows that the February diurnal wind velocity has a large amplitude 
(relative to the other periods investigated). The steepest velocity 
gradient occurs between 6:00 - 12:00. This is reflected in the o's 
computed around the pre-noon three hour average. 

From all the illustrated examples the January o's are the lowest. 
There are two relative maxima. Once in the morming at the 46 m level 


and one around noon-time at a lower level (%25 m). This is probably 


due to the small amplitude in the wind speed at the 46 m level (Figs. 13), ~ 


74 
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while at lower levels, the amplitude is larger and the variance increases. 


6. Stability characteristics of the forest and clearing 


An absolute measure of stability is non-existent. The Stability 
criteria are based, in most cases, On the available data. The most 
common stability criteria became the Richardson (Ri) number. 

The basic form of this parameter is the flux Richardson number 
defined as the ratio of buoyant production to stress production of 


turbulent kinetic energy 


ouw 93U/dz 


If the covariances are parameterized in terms of mean profiles 


- uw = vy, au/az 


- Ow = 36/az 


Vy 
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As the eddy viscosity and conductivity are unknown parameters, the 


23 


gradient Ri number is defined and used: 


The expression R_ = z_ 36/az 5 
Bo (2U/az) 


defines a local Ri number. The finite difference approximation, which 
is used for practical evaluation, is a bulk Ri number. ‘According to 
Deacon (1953), the time variation of the two is uniquely related at a 
fixed reference level. Therefore, the bulk Ri number can be used for 
profile classification purposes. 

Lettau (1957) used the geometric mean height for the height assign- 
ment of the difference quotient in wind and temperature. Dyer (1967) 
showed that the underestimation of the true gradient by the difference 
quotient, at the geometric mean height, is 2% for double levels and 
7.7% for fourfold intervals. 

To facilitate an approximate selection of the neutral cases to 
be used in computations of the micrometeorological parameters z,, d 


and uz, a bulk Richardson number was employed, suggested by Lettau 


(1957). 


(Ri)g = Z Ri(z5)/Z 2, 
j 3 
where 2, is the geometric mean height between two levels. (Ri), gives 
the overall convective stability and was found to be a useful criteria 


in this study. Lettau presents the following stability classification: 
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TABLE 4. STABILITY CLASSIFICATION OF BULK RICHARDSON GRADIENT NUMBER 
(After Lettau, 1957) 


= : -3 
Stability (Ri), , (10 37m) ss Stability (Ri), (1077 /m) 
"Class Class 
L3 = Extreme Lapse < than -15 Il = Weak Inversion 5 to 6 
12 = Moderate Lapse -15 to -10 I2 = Moderate Inversion 7 to 9 
Ll = Weak Lapse -9 to -5 13 = Strong Inversion 10 to 19 
N = Neutral -4 to 4 14 = Extreme Inversion > 19 


(Ri)p is only a gross indicator of stability above the forest 
canopy and in the two layers in the clearing. If a finer subdivision 
is made of the layer above the canopy, two adjacent sub-layers can have 
different stabilities. The plotted profiles suggested a separate 
classification of the following forest layers: 32 m - 36 m, 36m - 46 n, 
32 m- 46 m. For the clearing, the 1 m ~ 16 m and 16 m - 46 m layers 
were considered. The results are presented in Figs. 27,28,29. 
The different scales of the figures should be noted; the January 
values are an order of magnitude larger than the rest. 
A common feature in all the graphs is the diurnal variations of 
stability. During noontime it is negative, except for the 36 m - 46m ; 
layer in January. It seems that the (Ri) number, as represented in 
Fig. 19 allows a much better depiction of the avarage conditions in the 
whole layer than the use of only two measurement levels. Fig. 28 indicates 
that in spite of the diurnal variation in the stability parameter, all 
1 


the values are well within the -4 x 10° 2m -- 4x 10° 7a"! range, which 


was characterized by Lettau as neutral. 
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Fig. 29 represents the diurnal distribution of (Ri), number for 
the clearing. It is evident that a (Ri), number for the whole 1 m - 
46 m layer would be completely meaningless. The lower part is influenced 
by the strong insolation at the ground, the upper part being modified 


by the forest canopy. 


7. Determination of the micrometeorological parameters Za» d, uy 


For near neutral conditions it is hypothesized that the wind speed 
variation with height above the forest canopy is represented by 
u(z) = (v,/k) In [(z-d + 2,)/z5] 

i.e., the determination of the micrometeorological parameters zy »d 
and u, is necessary for profile fitting. The knowledge of Z, would also 
allow estimation of the energy dissipation in the planetary boundary 
layer (which is of synoptic importance). Information on the displace- 
ment height is necessary for aerodynamic heat flux computation, and 


uy gives the magnitude of the canopy stress. 


The factors that determine the roughness parameter are not yet 
fully understood. The interrelation of Z, and d with wind speed, above 


the forest, is still uncertain. The following represents a spectrum 
of opinions: 


Rauner (1960) - 2, increases and d decreases as wind increases 


Tajchman (1967) - 2g decreases and d remains constant as wind 
increases 


Allen (1968) - No dependence of z_ or d on wind speed 
Belt (1969) ° 
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Leonard & Federer (1973) - As Allen and Belt, when d is allowed to 
vary. When d is fixed, Z, shows a slight 
tendency to decline as wind speed increases, 
but only for west winds. 

Deacon (1953) - z, decreases with increasing wind speed. 

Kung (1963) - Vegetation cover rather than topography represents the 


most efficient roughness structure. 


7.1 Methods and results 


Graphical solution is frequently utilized to compute the micro- 
meteorological parameters Zo d and u,. The method is well justified 
for limited amount of data as an estimate. It involves quessing a 
plausible d value, plotting the result for each quess, and selecting 


the "best" fit. Many observed that the scatter of the z, values 


computed from different profiles could be quite large. Leonard and 
Federer (1973) give a range of 10-743 cm for a red pine forest (average 
tree height 11.6 m). Belt (1969) gives a range of 18-120 cm for a 
loblolly pine forest (average tree height 23.4 m). 

To obtain a representative average, it seemed necessary to use a 
large sample and therefore the graphical solution was impractical. 

The method which is being presently used is an adaptation of 
Stearns’ (1970) method. Surface roughness and displacement height are 


determined by an iterative technique, so that the error squares on wind 
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speed are minimum. E, is the error between the measured and theoretical 


wind speed, 


The two simultaneous conditions for a minimum are: 


3 2 
GZ = 22 (ZE;) 0 
Le] 
and 
C) 2 
GD = 3a CEs ) 0 


Given two initial guesses 2,41), Zo (2) and d(1), d(2) the rule of false 
position yields a new third approximation 
Z(1)62(2) - 2, (2)62(2) 


2,03) es ~ 6z(2) - Gza) 


where GZ(1) is evaluated using Z)(1) and d(1) as initial estimates of 
d. GZ(2) is evaluated with Z, (2) and d(l). If the tolerance check is 
not satisfactorily small, another approximation is made, replacing 
Z,(1) by 2,(2), GZ(1) by Gz(2), and 2,(2) by 2,(3). 

The final value of Z, together with d(1) is now used to determine 
GD(1). A second guess, d(2) is made and used and with 2,0) and Z, (2) 
for a new approximation to Z and subsequently GD(2). Then: 


d(1)GD(2) - d(2)GD(1) 
d(3) = GD(2) -  GD(1) 


ss 
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The iterative procedure was applied to the ensemble of 30-min averaged 
wind speed profiles in the 30 m - 46 m layer above the forest canopy, 
classified previously as neutral. The results of these computations 
are presented in Table 3. 

Z) and d estimates from measured profiles are used in literature 
to develop regression equations; especially equations which relate 
2 with vegetation height. Much caution should be exercised in using 
empirical regression equations. The constants in these expressions 
depend on such canopy characteristics as tree spacing, canopy shape and 
height. These equations are frequently used, because it is difficult 
to obtain reliable data for profile fitting. 

Stanhill (1969) gives the following relation between d and canopy 


height, h: 
log d = (.9793) log h - .1536 


This relation is in excellent agreement with values reported in the 
literature for vegetation canopies of lower height. 


Szeicz et al. (1969) found the following empirical relation between 


Zq and h: 
log (z,) = log (h) - .98 


Fig. 21 results from estimates of d and Zo using these formulae (for 
canopy height values of 29 m - 35 m), summary of several other studies, 
along with our results, in which Stearns' method was used. Stanhill's 
formula yields a d height some 4 - 10 m smaller than the values esti- 
mated from the profiles. The Szeicz formula yields a result which is 
consistent with the profile estimate for September but not for January 


or June. 


nae) 
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A test was run for possible correlation between Z and d with wind 
speed. Table 5 shows the results of a simple correlation-regression 


analysis applied to the ensemble of 384 cases in January, June, 


September. The z, is found to be an increasing function of wind speed, 


and d,a decreasing function of wind speed, which supports the findings 


of Rauner. Both parameters correlate best with the wind speed near the 


canopy. 
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SUMMARY 


This report presents the results of a comprehensive climatological 
and micrometeorological investigation of a dry tropical evergreen forest. 
The area under study is influenced by a cool, dry northeast monsoonal 
flow generally off the Asian continent during the period from November 
to February and by a warm, moist southwest monsoonal flow off the Andaman 
Sea and the Gulf of Thailand during the period from April to September. 
This latter period is characterized by two distinct rainfall surges in 
May and September separated by a period of lesser rainfall during June 
and July. The climatic regime at the forest site is somewhat similar 
to that at Bangkok, 190 km to the southwest. However, differences in 
elevation and significant differences in the heat and moisture balance 
in the forest yields a consistent cool bias of about 2.5°C in the 


forest and a more complicated interruption of the annual cycle of relative 


humidity compared with Bangkok. Both stations exhibit lowering temperatures 


and relative humidities during the cool northeast monsoon. The divergent 
flow around the Asian continental high pressure cell suppresses the 
precipitation and leads to an advective flux of moisture out of the 
Thailand region and relatively low precipitable water values. 

During the warm, moist southwest monsoon, the temperature, relative 
humidity, precipitable water and precipitation are all relatively high. 
Sakaerat, Bangkok and several other stations in southeast Asia show 
double maxima in precipitation during the season, the first surge 
during spring being supported largely by high amounts of evaporation 
but the second and larger surge during late summer and fall being sup- 
ported by a large advection of moisture into the region from maritime 


gources to the south and southwest. 
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The material presented in the report also depicts the seasonal 
variation of the microclimate of a dry evergreen tropical forest and 
a nearby clearing. Although the annual amplitude of ambient temperature 
in this part of the world is relatively small (approximately 6°C), the 
forest microclimate exhibits seasonal variability, induced primarily 
by the mode of heating and availability of moisture. The shielding 
effect of the canopy on the microclimate close to the ground in the 
“living space" is very much dependent on the prevailing cloudiness 
conditions which regulate the intensity of heating at the canopy top. 
It may be strong, amounting up to 4°C, or non-existent. 

Almost independent of the ambient wind speed the ventilation in the 
living space would be small, although the depth of penetration of wind 
flow below canopy varies. 

The relatively small clearing (500 m in diameter) was cut for 
the purpose of studying the effects of deforestation in this type of 
an environment. It was found that the surface temperature difference 
between the clearing and below the canopy during the day can reach 
up to 46°C. During the night, the clearing loses heat more effectively 
than the forest floor, resulting in temperatures which are closer 
to each other. 

The forest has a significant effect on the wind field reduction, 
in particular, during weak flow. Below the canopy top, the diurnal 
variation is almost non-existent. In the clearing, the winds close 
to the ground do not exhibit a strong diurnal variation, except during 


January and December. 
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For an assessment of the mean stability characteristics of the forest 
site, values of the bulk gradient Richardson number were computed. 

A frequency distribution of the number of occurrences in each stability 
class was obtained. Nearly all the cases sampled for the above forest 
canopy during the moist, southwest monsoon period fell into the near- 
neutral classification. For the same layer in January, during the cool 
dry northeast monsoon, about 41% of the cases exhibited near-neutral 
conditions with slightly less than 10% falling in various of the 
unstable categories and slightly less than half falling in various of 
the stable categories. For the clearing, the lower layer exhibits a 
wider variation of stability than does the upper layer. Slightly less 
than 80% of the cases in June and about 70% of the cases in September 
show near-neutral stability in the upper layer of the clearing; evidence 
of the modifying influence of the surrounding forest. 

For near neutral conditions in both monsoon periods, the winds 
in the microlayer above the forest canopy show a reasonable fit to a 
displaced logarithmic profile. The datum displacement height has an 
average value of about 27.6 meters but varies over a range of about 
2(%) meters from the northeast to southwest monsoon periods. 

The roughness parameter for the canopy exhibits a strong seasonal 
variation with a mean value of 0.83 m for the January sample, 5.55 m 
for the dene sample and 3.61 m for the September sample. The canopy 
stress has a mean value of 1.3 dynes-cm~2 for January and the much 
larger mean values of 21.6 dynes-cm72 for June and 13.7 dynes-cm~2 


for September. 
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A test was made for possible correlation between z, and d with 


oO 
wind speed at the top of the tower (46 m). The roughness parameter 

was found to be an increasing function of wind speed and the datum 
displacement height a decreasing function of wind speed. Both parameters 
correlate best with the wind speed near the canopy (i.e. u(32)), but 


the correlation between d and u(46) is also significant at the 5% level 


of significance. 


APPENDIX A 


Present "D" Tapes Status 


Present "D" Tapes Status 


In this section a report is provided on the present status of 
the TREND "Dp" data tapes at the University of Maryland. 

All the tapes at our disposal are copies of the original tapes 
taken at the Thailand forest site. The entire set is comprised of 
two types of tapes; 1) digital tapes ("D'"' tapes) - containing infor- 
mation about temperature and wind velocity, measured every 10 seconds 
2) analog tapes ("A'' tapes) - containing information about radiation, 
humidity and wind direction, measured every 30 seconds. Part of the 
"D" tapes were reduced and packed in binary. A "Users Manual" explaining 
how to use the repacked binary tapes is in preparation. 

In Table Al the numerical code of the packed tapes and the cor- 


responding final product are presented. 


Table Al 

No. Packed "D" Tapes 

240 237, 238, 240, 241, 243 
275 273, 274, 275, 277, 278 
287 285, 286, 287, 288, 290 
293 291, 292, 293, 295, 296 
301 298, 299, 301, 302, 304 
308 305, 306, 308, 309, 310 
316 312, 314, 316, 317, 318 
331 328, 330, 331, 333, 334 
338 335, 337, 338, 340, 341 
354 351, 352, 354, 355, 356 
363 363, 365, 366, 367, 369 
374 

380 377, 378, 380, 381, 383 


395 393, 394, 395, 397, 398 
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In Table A2 the conditions and dates covered by each of the reduced 
tapes are indicated. It is evident from Table A2 that the useful in- 
formation retrieved from the raw data is not continuous, but rather 


in stretches of approximately two weeks length. 


TABLE A2 


"Dp" Data Tapes 


Dates Remarks 
240 Sept. 25 - Oct. 8, 1969 99% good -ecans 
275 Dec. 1 - Dec. 12, 1969 good 
287 Dec. 26. - Jan. 8, 1970 good 
283 Jan. 8 - Jan. 19, 1970 very good 
301 Jan. 19 - Feb. 15, 1970 good (gap in data for Jan. 
21 - Feb. 5, 1970) 
308 Feb. 15 - Feb. 28, 1970 very good 
316 Feb. 28 - March 14, 1970 good 
331 ; April 6 - April 20, 1970 very good 
338 April 20 - May 5, 1970 good 
354 May 24 - June 6, 1970 80% good 
363 June 19 - July 2, 1970 99% good 
374 July 2 - July 24, 1970 very good 
380 July 24 - August 11, 1970 ? 
395 August 29 - Sept. 13, 1970 75% good 


Additional "D" tapes are at our disposal. These tapes were not yet 


reduced because of the cost involved. Approximately, 105 tapes were 


not reduced. 
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APPENDIX B 


Description of the condensed 
data tape 
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From the overwhelming data base at our disposal, a "working" data 
set was prepared. This unit consists of half hourly averaged data obtained 
from the densely packed binary tapes of the ten second data. These half 
hourly averages were obtained for various periods of the monsoon cycle, 
for about two weeks each. Since the reduction process was time consuming 
and expensive, this useful set of data was stored on a separate tape. 

From the half hourly data base described previously, two week averages 
were obtained. These are also stored on the same tape. 

In this section, information about the tape content, and format 


will be presented. 


Instructions for condensed tape users. 


This 9-track tape (intemal tape Number P10629) contains 170 files, 
each of which contain wind speed and temperature data. The tape is 
unlabled and contains fieldata code character representation with a 
storage density of 6250 bpi. 

Each file is 2 blocks long with the first file starting at the load 
point of the tape. 

To access any file the following procedure is suggested for use 
on Univac 1100 operating system: 


Assign the tape: 
@ASG,TJ TAPE,U9S,p10629 


Assign an alternate mass storage file, NORMAL CONVENTIONS: 
@ASG FILE 


Move the tape to the desired file; (refer to lisbd: 
@MOVE tape., File # 


Copy the tape file to the mass storage file: 
@COPY Tape.,File 
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After the copy is produced the tape will be positioned at the beginning 
of the next file (after the EOF mark). Care must be taken to ensure 
you don't @MOVE past the end of the tape. 

Now the mass storage file can be used as a standard Univac data 
file. 

The file will consist of 48 records each of which is a hourly 


average (240 cards total); i.e., 5 cards per record: 


1 - Pate and time 
2,3,4 - 21 temperature levels 
5 - 8 velocity levels 


First will come the clearing tower (CT) data followed by the forest 
tower (FT) data. 
Card Formats: 
In general the card format will be as follows 
14,4X,13,1X,212/10(F6.2,1X) /10(F6.2,1X) /F6.2/8(F5.2,1X) 
when noted by ('*') the card format will be 
3(3(12,1X),1X)/10F(6.2,1X)/10(F6.2,1X) /F6.2/8(F5.2,1X) 
(the only difference being in the date format). 
The following list describes the content of each file stored on 


the tape. All the data stored are for 1970 except December (1969). 


a) 


File # 


The approximately two week averaged data 


The 5 hourly averaged data for single days 


156 
157-162 
163 
164-169 


Contents 


Dec. Avg. 
Jan. Avg. 
Feb. Avg. 
Apr. Avg. 
Jun. Avg. 
Jul. Avg. 
Sep. Avg. 
Dec. Avg. 
Jan. Avg. 
Feb. Avg. 
Apr. Avg. 
Jun. Avg. 
Jul. Avg. 
Sep. Avg. 


(Julian date) 


336-338 
341-345 
336-338 
341-345 
4-18 
4-18 
47-58 
47-58 
97-98 
100-109 
97-98 
100-109 
171-179 
181-182 
171-179 
181-182 
183-184 
186-190 
200-204 
183~184 
186-190 
200-204 
224 
224 
247 
249-254 
247 
249-254 


(4-8) 
(4-8*) 


Information about missing data 


Several profiles during yaetote periods of time were missing. 
In order to obtain a complete set of data, consisting of 48 half hourly 
profiles for each day, limited editing was exercised. The following 
procedure was followed. Each missing profile was replaced by a profile 
from a complete day, when there was a good agreement between the data 
preceeding and following the missing profile. 


The following list presents a full documentation of the data that 


were supplemented. 


The replaced data in the files are as follows: 


Missing Replaced by Where 
Day Year Time Day - Year Time CT&¥FT 
338 1969 10:00 338 1969 10:30 FT&CT 
342 1969 20:00 342 1969 20:30 FT&CT 
345 1969 8:00 345 1969 7:30 FT&CT 
345 1969 8:30 345 1969 9:00 FT&CT 
7 1970 15:00-16:30 8 1970 15 :00-16:30 FT&CT 

11 1970 3:00- 6:30 ll 1970 2:30,2:30, 

2:30,2:30, 

7:00,7:00, 
7:00,7:00 FT&CT 
175 1970 0:00 1970 cT 
177 1970 7:00 177 1970 6:00 CT 
178 1970 11:30 178 = 1970 11:00 CT&FT 
181 1970 0:00 179 1970 23:30 CT&FT 
183 9:30 180 1970 9:30 FT 
186 0:00- 1:30 185 1970 0:00- 1:30 FT 
187 18:00-23:30 186 1970 18:00-23:30 FT 
188 0:00-13:00 187 1970 0:00-13:00 FT 
200 13:30-16:00 201 1970 13:30-16:00 FT 
202 12:00-18:30 201 1970 12 :00-18:30 FT 
204 8:30 208 1970 8:30 FT 
247 14:30-16:00 253 1970 14:30-16:00 cT 
250 23:00 246 1970 23:00 CT 
252 20: 00-21: 30 254 1970 20:00-21:30 CT 
247 14: 30~-16:00 251 1970 14:30-16:00 FT 
250 23:00 252 1970 23:00 FT 
253 20:00-21:30 254 1970 20:00-21:30 FT 
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Table 2. Two-week averaged wind velocities at selected levels 
(V1 ~- the 46 m level Vy - the 32 m level) 


Period 
No. of Convergent 
Cases 


January &-10,1970 
40 cases 


June 20-26,1970(A) 
201 cases 


June 20-26 ,1970(B) 
72 cases 


Sept, 9-10, 1970 
67 cases 


July 19-23, 1970 
108 cases 


August 11-13, 1970 
68 cases 


Averag2 (weizhted 
by No. of 
cases) 


4.25 
(.89} 


5.41 
[1.9] 


4.69 


d(m) 


{o, (a) ] 


29.53 
[ 0.05) 


27.18 
{ 0.69} 


27.16 
{ 0.72) 


28.09 
{ 0.75) 


27.49 
[.47] 


27.67 
[.85) 


; 27.58 


1.20 
[.43] 


1.41 
{.48] 


1.21 


t(u,) (dynes grey 


1.31 


24.2 


19.1 


13.7 


17.28 


23.86 


18.97 


Table 3. Average values and standard deviations of roughness 
parameter, datum displacement height, friction velocity 


and canopy stress for the Sakaerat Forest. 


{Xd + V = 4] 
; 3so10q pueT}eyL ay2 aaoqge 1zaquaqdas 
pue sun¢ ‘Azenuer uy satTytjord puytm pa8ezaae aynuyu-o¢ 


JO S9SBD HEE oj} SOTIST~IeIS UOTSseI9e1 pur UoTIeT e110) "¢ aTqey 
ay2 38 abcseuaes . 

yzord = 9S OLL*0 9£000°0~ B86°L 82 Tet Leo°a e276" (ze) [ (wa) PUT 

820°0 4682 9$9°O £T000°0~ £66°¢ 6€T 627 LE0°O €26°L (99)" ( (a>) P) BT 

70S°O nce 702°0 ™900°0 808"? el Tet 80L°0 896°S (7¢)9 ( (wo) °z]ut 

£95°0 §=6 €2Z,—s« L09°0 g0f£00°0 $99°9 6£T 677 B0L°0 896°S (oyyn (mo) °z] UT 

679 = 89S ML L"0 LT0°T= $762 BL 18t 20T 192 (ze)A (=>)p 

Ole C62 659°0 €87°0- 8962 6€T 627 20T 1942 (99)" (w2)P 

B'79T C92 = 6E9O esi*T pce 82 TST 912 999 (20)8 (m>)°2 

S°CLT 002 985°0 006°0 €°6L 6£T 629 912 999 (9y)8 Coe 

ress ial Tx/a) Cay [ 995-83) (Tes TAT 293 an + SLINN 
aaBwt sy JO Ofer (UOT Ze{I2109) io = 7 cy are 
2012g pawpurys i Y @ v oe x 4 x x k 


en = Se etn er seen. 


*- 


Tat meine 


APPENDIX C 


Listing of selected data 


‘ss 


The data presented in this section represent the following two week 


averages: 


December (1969), January, April, June, 


July, September (1970) 


The data are presented first for the forest and then for the clearing 


in the following format: 


a) 


b) 


c) 


average wind speeds at eight levels* 


average temperatures at fifteen levels 
above the surface* 


average temperatures at six levels - 
surface and subsurface 


*For details about the location of the various levels see Fig. 3 
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